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Abstract 

Gribov theory is applied to investigate the shadowing effects in the structure functions 

^ , of nuclei. In this approach these effects are related to the process of diffractive dissociation 
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1. Introduction 

Deep inelastic scattering (DIS) on nuclei gives important information on distributions 
of quarks and gluons in nuclei. The region of small Bjorken x is especially interesting 
because partonic clouds of different nucleons overlap as x ^ and shadowing effects 
become important. There are experimental results in this region, which show that there 
are strong deviations from an A'^ behavior in the structure functions [1] . Several theoretical 
models have been proposed to understand these data [1]. The most general approach is 
based on the Gribov method [2]. It relates partonic and hadronic descriptions of small 
X phenomena in interactions of real or virtual photons with nuclei. In this approach the 
shadowing effects can be expressed in terms of the cross-sections for diffraction dissociation 
of a photon on a nucleon (Fig. 1). This process has been studied recently in DIS at HERA 
[3]. The detailed x, Q^ and M^ (M is the invariant mass of the diffractively produced 
system) dependencies observed in these experiments [3] have been well described in the 
theoretical model of ref. [4] which is based on Regge factorizations and uses as an input 
available information on diffractive production in hadronic interactions. Here we will apply 
the same model to calculate the structure functions of nuclei in the small x-region. The 
use of the model, which describes well the diffraction dissociation of virtual photons on a 
nucleon target, leads to a strong reduction of the theoretical uncertainty in calculations of 
the structure functions of nuclei in comparison with previous calculations [1, 5]. 

2. The model 

In the Gribov approach the forward scattering amplitude of a photon with virtuality 
Q^ on a nuclear target can be written as the sum of the diagrams shown in Fig. 2. The 
diagram of Fig. 2a corresponds to the sum of interactions with individual nucleons and is 
propotional to A^. The second diagram (2b) contains a double scattering with two target 
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nucleons. It gives a negative contribution to the total cross-section, proportional to A^'^ 
(for large A) and describes the shadowing effects for sea quarks. According to reggeon 
diagram technique [6] and Abramovsky, Gribov, Kancheli (AGK) cutting rules [7], the 
contribution of the diagram of Fig. 2b to the total 'j*-A cross-section is related to the 
diffractive production of hadrons by a virtual photon as follows : 



r r drr'^ 



FA{tmin) 1 (1) 



t=0 



where Ta(6) is the nuclear profile function, pA is the nucleon density (Ta(6) = /_^ dZ pA{h, Z), 
j(fhTA{h) = 1) and 

FAitrmu) = J dZ e-"^^^pA{b, Z)/TA{h) , trmn = -qI = -m% X^ f ^,^ J 

Note that FA{tmin) is equal to unity as x ^ and decreases fast as x increases to Xcr ~ 
— -^p— , due to a lack of coherence for x > Xcr- 

Thus the second order rescattering term can be calculated if the differential cross- 
section for diffractive production by a virtual photon is known. 

Higher order rescatterings are model dependent, but in the region of x>10~"^, where 
experimental data exist, their contribution is rather small (~ several %) for existing nuclei. 
We use the following unitary expression for the total 7*-A cross-section 



where 



f{x,Q^) =4n f dM^ 



dai^,p 



FA{tmin)/cf-i*N 
t=0 



dW^dt 

This expression is valid in the generalized Schwimmer model [8, 9] and corresponds to the 
summation of fan diagrams with triple Pomeron interaction. We have checked that the 
results for eikonal-type summation of higher order rescatterings are very similar to the one 
obtained with (2). 



Thus we obtain for the ratio Ra = F2A/F2N, in the region of smaU x 

^- - /■,.,^^^4^l(^)^ . (3) 



F^K J l + (A-l)f{x.Q-2)TA{b) ' 
The deviation of this ratio from A^ = A j d'^bTA{b) is due to the second term in the 
denominator of the integrand in eq. (3). Thus, knowing the differential cross-section for 
diffraction dissociation on a nucleon and the structure function of a nucleon (cr-y.Ar), one 
can predict the A (and x, Q'^) dependence of structure functions of nuclei. Eq. (3) can 
only be used in the region x < 10~^ where the sea quarks component dominates, because 
at larger values of x shadowing of valence quarks (which in general is not described by 
eq. (3)) becomes important [10, 11]. The effects which lead to antishadowing are also 
important in the region of x ~ 0.1. 

In refs [4] we described the diffractive contribution to DIS in terms of Pomeron ex- 
change 

F^ix^Q^xp^t) = ^-^^P^x'-'-^^'^FpiP^Q^t) (4) 

Idtt 

where gppit) is the Pomeron-proton coupling [gppit] = (7^(0) exp(Ct) with ((7^(0))^ = 

23 mb and C = 2.2 GeV'^), ap{t) is the Pomeron trajectory (ap(0) = 1.13, a'p{0) = 0.25 

GeV~^) and Fp{/3,Q'^,t) is the Pomeron structure function. The variable /3 = j^i.Qi = 

-^- plays the same role for the Pomeron as the Bjorken variable x for the proton. At large 

Q^, Fp can be expressed in terms of the quark distribution in the Pomeron 

Fp(AQ^t) = 5^e?/3 [qf{(3.Q^t) + q[{(3,Q^t)] . (5) 

i 

In refs. [4] we determined Fp{(3j Q^, t) using Regge- factorization for small values of (3 and 
a plausible assumption on the /3 -^ 1 behavior. This function was then used as an initial 
condition for QCD evolution of partons in the Pomeron. The results of the QCD-evolution 
crucially depend on the form of the gluon distribution in the Pomeron. Experimental 
results for F^ can be understood only if the distribution of gluons in the Pomeron is 
rather hard and the gluons carry the main part of the Pomeron momentum [4, 12-14]. 



The model described above gives a reasonable description of the process of diffractive 
production in DIS and can be used to compute the function /(x, Q^), which determine the 
shadowing of nuclear structure function via (3). This function can be written in terms of 
the ratio Fp/ F2N '■ 

fix, Q) = jj-^ fep(O)) [-) FMi3;Q2-fA{t^^n) (6) 

where the integral limits are x/xqp with xqp = 0.1 and Q'^/{m'^ + Q^). For the nucleon 
structure function F2N we use the results of ref. [15] and for Fp{(3,Q'^) the expression 
and parameters of ref. [4b]. In the numerical calculations we use a standard Woods- 
Saxon profile for A > 20 and a gaussian for A < 20 with a rms radius given by Ra = 
0.82 A^'^ + 0.58 fm [16]. For deuteron we have also used a gaussian with an artificially 
high radius in order to reproduce the small amount of shadowing present in this case. The 
results of our calculations are shown in Figs. 3-7. 

3. Numerical results 

Comparison of our predictions for the ratio ^F2a/F'P with experimental data of 
NMC [17] is shown in Fig. 3 and for ratios of different nuclei in Fig. 4. New data for the 
ratio Fi-/F^ [18] are also shown in Fig. 4. It is important to note that experimental 
points in Figs. 3, 4 for different x correspond to different values of < Q^ > [17] [18]. 
This correlation has been taken into account in our calculations. The agreement between 
theoretical predictions and experimental data is rather good taking into account that there 
are no free parameter in our calculations. However, one should notice that for large nuclei 
(Sn), the onset of the shadowing seems to move to smaller values of x as compared to 
light ones. This trend is also present in the Pb data of ref. [19] and is not reproduced in 
our model (despite the decrease of Xcrit with increasing A) . This probably indicates that 
antishadowing (not present in our model) extends to smaller x values for heavier nuclei. 

Our predictions for the ratio ^ jr^ in the region of very small x are shown in Figs. 5 
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for fixed values of Q^ . They can be confronted to experiment if nuclei at HERA will be 
available. Note that our results are more reliable for small values oi x [x < 10~^) where 
the effects of both valence quark shadowing and antishadowing are negligeable. The curves 
for shadowing effects in the gluon distribution of nuclei 4- ^ N(n-i\ ^^^ shown in Figs. 6. 
They look similar to the shadowing for the quark case. However the absolute magnitude 
of the shadowing is smaller in the gluon case contrary to expectations of some theoretical 
models [1] but in agreement with [20]. (Note that these results are sensitive to the gluonic 
distribution in the Pomeron, which is poorly known at present). These predictions can be 
tested in experimental studies of J/ifj and T-production on nuclear targets at RHIC and 
LHC. 

Finally we want to discuss in more detail the Q^-dependence of the shadowing. Recent 
NMC data [18] for the ratio of F2 "/-F^l" ^^^ shown in Figs. 7 as functions of Q^ for fixed 
values of x in the small x region. The theoretical curves have a weak dependence on Q^ 
and are in a reasonable agreement with experiment, although the Q^ dependence seems 
stronger in the data especially in the region of small Q^. At larger values of x the data are 
practically Q^-independent. These properties should be checked in future experiments. 

4. Conclusions 

In conclusion, a model based on the Gribov-Glauber theory of nuclear shadowing 
and the properties of diffraction in DIS observed at HERA, leads to a fair description of 
experimental data on structure functions of nuclei in the small x region. Predictions of 
shadowing effects for quark and gluon distributions are given. They can be tested in future 
experiments at HERA and in hadronic colliders. 
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Figure Captions 

Fig. 1 : Diffractive dissociation of a virtual plioton. Tlie sliaded area represents tlie 
exctiange of a Pomeron. 

Fig. 2 : Tlie first two terms (single and double scattering) of ttie multiple scattering 
series for the total 7*A^ cross-section in the Gribov- Glauber approach. 

Fig. 3 : The ratios {2/A)F^/F^ computed from eq. (3) for different values of x. The 
experimental points are from ref. [17] . The values of Q^ are different for different 
x- values [17]. 

Fig. 4 : The ratios {Ai/A2)F2 '^ jF^ ^ computed from eq. (3) for different values of x. 
The experimental points are from refs. [17] and [18]. The values of Q^ are 
different for different x values [17, 18]. 

Fig. 5 : The ratios {Ij A)F2 IF2 computed from eq. (3) for different values of x in the 
small X region, at fixed values of Q^ . 

Fig. 6 : The ratios {Xj A'i)g j g^ of gluon distribution functions computed for different 
values of x in the low x region, at fixed values of Q^ . 

Fig. 7 : The ratio (12/119)F2^"/F2^ computed from eq. (3) for different values of Q^, 
at two fixed values of x. The data points are from ref. [18]. 



1^2 /^2 



4 F,^ 



2^Fl_ 
12 F^ 



2_f£^ 

40 F^ 



1.1 

1.05 

1 

0.95 
0.9 

0.85 
0.8 

1.1 

1.05 

1 

0.95 

0.9 
0.85 

0.8 

1.1 

1.05 

1 

0.95 

0.9 
0.85 

0.8 



-1 1 1 1 1 r- 



-1 1 1 1 1 1 1 [— 



He/D 




H 1 1 1 1 — h 



C/D 



H 1 1 1 1 — h 



Ca/D 



H 1 1 1 1 1 — h 




H 1 1 1 1 1 — h 




0.01 



0.1 



Fig. 3 



^2 Z / " ^ 



±11. 

12 F^' 



6 Ff" 
40 i^2^^ 



1.1 

1.05 

1 

0.95 
0.9 

0.85 
0.8 

1.1 

1.05 

1 

0.95 

0.9 
0.85 

0.8 



12 ^2 
40 F^ 



1.1 
1.05 

1 - 

0.95 - 

0.9 - 

0.85 - 

0.8 - 



12 ^2 

119 F^ 



1.1 

1.05 h 

1 
0.95 h 

0.9 
0.85 

0.8 h 
0.75 

0.7 - 

0.01 



-1 1 1 1 1 r- 



C/Li 



J— i-r^— ^— r^-^ 



T I 



1±_^^^^^ 



Ca/Li 



Ca/C 



Sn/C 



H 1 1 1 1 h 



i I I !_!_ 



H 1 1 1 1 h 




H 1 1 1 1 h 



0.1 



Fig. 4 



(5' - h.OGeV^ 



2 o.s 



AFi' 




H 1 1 — I I I I I I 1 1 1 I I I I I I 1 1 1 I I I I I 

Q^ = 25.0GeV^ 



H 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I 

Q^ = lOO.OGey^ 



AFi^ 




Fig. 5 



= S.OGeV^ 



AgN 




H 1 1 — I I I I I I 1 1 1 I I I I I I 1 1 1 I I I I I 

Q^ = 25.0GeV^ 



Q^ = lO.OGey' 



H 1 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I 



Q^ = lOO.OGeF^ 



AgN 




Fig. 6 



12 TpSn I TpC 
119^2 /^2 



0.98 


- 


1 1 




1 


1 1 1 1 1 


0.96 


- 








- 


0.94 


- 










- 


0.92 

12 i^9^" 
^t£ 2^ 0.9 

119 Ff 

0.88 


- 


J J 


J 


I 


-- 


- 


0.86 


- 








- 


0.84 


- 


-- 






- 


0.82 


- 








X = 0.0125 


0.8 












0.98 


- 










- 


0.96 


- 










- 


0.94 
0.92 

12 i^9^" 
^t£ 2^ 0.9 

119 fC' 

"^ 0.88 


- 






J 


■■ 


■■ 


I . 


■ - 




^ 






- 


0.86 


- 








- 


0.84 


- 








- 


0.82 
0.8 


- 








a; = 0.0175 



10 



Q' 



Fig. 7 




A 




Fig. 1 





(a) 



(b) 



Fig. 2 



